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Exact static solutions to a translationally invariant discrete 4 model 
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For a discrete, translationally-invariant (f> 4 model introduced by Barashenkov et al. [Phys. Rev. 
E 72, 35602R (2005)], we provide the momentum conservation law and demonstrate how the first 
integral of the static version of the discrete model can be constructed from a Jacobi elliptic function 
(JEF) solution. The first integral can be written in the form of a nonlinear map from which any 
static solution supported by the model can be constructed. A set of JEF solutions, including the 
staggered ones, is derived. We also report on the stability analysis for the static bounded solutions 
and exemplify the dynamical behavior of the unstable solutions. This work provides a road-map, 
through this illustrative example, on how to fully analyze translationally-invariant models in terms 
of their static problem, its first integral, their full-set of static solutions and associated conservation 
laws. 

PACS numbers: 05.45.-a, 05.45.Yv, 63.20.-e 



I. INTRODUCTION AND SETUP 

In recent years, the discrete nonlinear models have 
played a very important role in several physics applica- 
tions P, E| ■ The question of mobility of solitonic exci- 
tations in discrete media is a key issue in many physical 
contexts; for example the mobility of dislocations, a kind 
of topological solitons, is of importance in the physics of 
plastic deformation of metals and other crystaline bodies 
3] . Similar questions arise in optics for light pulses mov- 
ing in optical waveguides or in photorefractive crystal 
lattices (see e.g., [|| for a relevant recent discussion) and 
in atomic physics for Bose-Einstein condensates moving 
through optical lattice potentials (see e.g., [|[ for a re- 
cent review). These issues may prove critical in aspects 
related to the guidance and manipulation of coherent, 
nonlinear wavepackets in solid-state, atomic and optical 
physics applications. 

In the above settings, the translationally invariant (TI) 
discrete models [|| have received considerable attention 
since they admit static solutions that can be placed any- 
where with respect to the lattice. For the Hamiltonian TI 
lattices 0,11], this can be interpreted as the absence of the 
Peierls-Nabarro potential For the non-Hamiltonian 
lattices, the height of the Peierls-Nabarro barrier is path- 
dependent but there exists a continuous path along which 
the work required for a quasi-static shift of the solution 
along the lattice is zero [9(. In general, one can state that 
coherent structures in the TI models are not trapped by 
the lattice and they can be accelerated by even a weak 
external field. This particular property makes the TI 
discrete models potentially interesting for physical appli- 
cations and one such physically meaningful model has 
been recently reported [l(| ■ 

It has been demonstrated that some of the TI lattices 
support traveling solutions, but only for selected veloc- 
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completely integrable lattices 0, [H, [3] where the prop- 
agation velocity can change continuously. Translationally 
invariant discrete models can also be defined as models 
whose static version is exactly solvable; thus, the TI dis- 
crete models are directly connected to the theory of inte- 
grable maps [lH ]. Starting from the pioneering works of 
0, 0] , the TI discrete models have been receiving atten- 
tion from several groups. General classes of such models 
have been constructed and investigated for the Klein- 
Gordon @, i, M, [U GJ, Ell Hand for the nonlinear 
Schrodinger equations (NLSE) GJ, M, EH El SI • 

For a wide class of Klein-Gordon TI models it is known 
that they conserve momentum defined in a standard form 
[f| but are non-Hamiltonian. That is to say, there are TI 
models that conserve the standard momentum but not 
the standard energy. Similarly, there are models that 
conserve energy but do not conserve the momentum; see 
{16], for a detailed discussion. However, this raises the 
question: is there always some form of a conservation 
law (either as a standard momentum or a standard en- 
ergy, or as a generalized version thereof) present in the 
TI models? In this paper, we consider an example of a 
model that conserves neither the standard energy, nor the 
standard momentum. We show in this case that a conser- 
vation of a non-standard momentum is present. The ex- 
istence of a conservation law even in such a non-standard 
case, leads us to believe that there should always be at 
least one conservation law in the TI models. 

It is also known that some of the discrete TI models, 
both of the Klein-Gordon and of the NLSE type, support 
exact static/stationary and sometimes even moving (with 
a discrete set of velocities) solutions in the form of the 
Jacobi elliptic functions (JEF) [1, ©, GJ, EH El but other 
TI lattices do not support JEF solutions, e.g., the lattice 
reported in [3]. 

On the other hand, for any TI lattice, the static prob- 
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lem can be reduced to a first-order difference equation 
which has been called the discretized first integral (DFI) 
0, E3| because, in the continuum limit, the DFI reduces 
to the first integral of the static field equation. The DFI 
is a nonlinear map from which, for an admissible initial 
value, one can construct a static solution iteratively. The 
DFI approach can be used to construct the TI models 
@i E3 , but here we address the inverse problem of find- 
ing the DFI for a given TI model. We demonstrate how 
the DFI can be found from a known JEF solution to a 
TI model. 

More specifically, in the present study we consider the 
following discrete model first suggested by Barashenkov 
et al. [3, 
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ft2 



2(j) n + (j)n+l) + A</>„ - \(f) n -i(f> n (f) n+1 ,(l) 



where A = ±1 and overdot means derivative with respect 
to time. In the continuum limit (h — > 0) the model re- 
duces to the <fi 4 field theory model 



+ Ai 



(2) 



The static version of Eq. @ (with omitted inertia term 
4>u) possesses the first integral 
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l ) 2 + K = 0, 



(3) 



where K is the integration constant. An exact static 
kink solution to the TI lattice of Eq. (JT|) and a two- 
point map that generates the kink solution were derived 
by Barashenkov et al. [Hj]. 

There are numerous reasons for the choice of Eq. (Q]) 
as our illustrative example for this study, intending to 
highlight the range and the power of the methods avail- 
able to tackle TI models. Firstly, the model does not 
possess an obvious conservation law. Secondly, it does 
not have an obvious discrete first integral. Thirdly, its 
general solutions are not available and neither their para- 
metric space, nor their stability have been explored. In 
the present work, we develop a number of relevant results 
for the lattice of Eq. ([I]). In particular, in Sec. HH we give 
(i) the corresponding conservation law; (ii) the two-point 
map (or DFI) from which any static solution of Eq. ((T|) 
can be constructed; (iii) a set of static JEF solutions. We 
believe that the features (i) and (ii) should be present in 
any TI model, and use this non-standard model to show- 
case how to develop them. Then, in Sec. IIIH we present 
the results of numerical studies including the stability 
analysis of the static solutions and, for h 2 < 2, we show 
several examples of the dynamics of unstable solutions. 
Section [TVl concludes the paper. It is worth pointing out 
that recently two of us [§| have similarly studied a model 
(henceforth called the -standard- momentum conserving 
or MC model) 
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2</>„ + (f> n+ i)+\(j>n—^r4'n[ ( i > n-l+ ( l ) 



We shall compare and contrast the various results ob- 
tained from model of Eq. (Q]) with those of the MC model 
to get a better insight into the various issues involved. 



II. ANALYTICAL RESULTS 

A. Momentum conservation law 

By seeking expressions that would symmetrize the 
form of the nonlinearity, we have found that the discrete 
TI model of Eq. {T]) conserves the momentum defined as 



P = ^2 ^n{<f>n+2 - 4>n-2), 



(•5) 



which is different from the standard definition of momen- 
tum used, e.g., in [a. Ila]. The latter momentum, namely: 



>n\9n+l ~ (Pn-1), 



(G) 



is conserved e.g. by the MC model of Eq. (j4|). The 
proof in both cases is through straightforward time- 
differentiation and application of a telescopic summation. 



B. Reduction to the first-order difference equation 

It is possible to construct the first integral of the static 
version of Eq. ([1]) , which can be presented in the form of 
the nonlinear map, 
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Z(j> n ± f{4>n) 



K-X<j>l 



with 



A = A/i 2 , 



K = 1 



2K 



z = 



(2 - A) 2 - KA 2 
2(2- A) 



X 



KA 2 + (2 - A) 2 (l - Z 2 ) 
A(2- A) ' 



(7) 



(8) 



(9) 



(10) 



(11) 



(4) 



Apart from the model parameters, h 2 and A, the map 
defined by Eqs. ([7]) to (fTTj) contains the integration con- 
stant K. The role of the second integration constant, 
for the second-order difference static equation ([TJ) , is ef- 
fectively played by the initial value of the map, 4>o- I n 
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Appendix I we show how the two-point map can be con- 
structed from the known JEF solutions. Due to the sym- 
metry of Eq. JTJ), m the left-hand side of Eq. |[7J), one 
can substitute 4> n ±i with 4> n ^i- 

For any pair of admissible values, K and (f>o, the map of 
Eq. (O generates a static solution to Eq. (UJ. Since the 
same map is used for calculation of the forth and back 
points of the map, iterations starting from an admissi- 
ble value (j)o cannot give an inadmissible one and thus, 
the static solution will necessarily be constructed for the 
entire chain. From the two different roots of Eq. ([7]) 
for (j) n +i (or for </>„_i), one should take the one different 
from 4> n -\ (or from (f> n+ i). This rule gives the sufficient 
condition for moving along the same solution branch and 
it guarantees that, for the chosen (f> n +i (or for </>„_i), the 
three point problem Eq. |T]) will be satisfied for </>„. 

In the continuum limit, the first-order difference equa- 
tion ([7]) reduces to Eq. ([3]). In other words, the two-point 
map, as given by Eq. jJJ, is a DFI of the static (j) 4 field. 
This observation is in line with our earlier work [l7| 
where the relation between TI discrete models and DFI 
of the static field equation was established. 

It is interesting to note that for the MC model the 
nonlinear map has the form Q 



(12) 



(2 - A)0 ?1+1 ± y/D 



where 



D = 2A(1 



(2- 



J n+1) 



A^ +1 ) 



2C*(A0 



n+1 - 2) 



(13) 



and C = Ch 2 with C being the integration constant. 

It is also worth pointing out that while the map for the 
MC case given by Eqs. (fT2j) and (fT3|) . is quadratic in both 
4> n and </> n +i, the map in the present case as given by Eq. 
([Jj is, generally speaking, quartic in <p n and quadratic in 
4> n +i, the only exception being when A = 2 when the 
map is quadratic in both <j) n and (this is discussed 

in detail in Sec. Ill C 4j) . Thus, while the map in the MC 
case belongs to the family of twelve-parameter two-point 
map presented by Quispel et al. [lj| , the map in our case 
as given by Eq. (UJ does not, in general, belong to the 
relevant family of maps. To be more precise, while the 
static version of Eq. jT]) is of the Quispel et al. [UJ form 
as given by 



h{4> n )-h 



M<M - M 



9n-l 



(14) 



with 



/i(0„) = (2 - K)4> n , f 2 {cp n ) = 1 , h{4>l) = K4> n ,(15) 

one is unable to write the corresponding invariant map 
in terms of the twelve parameters presented in [TfJ. 

It should however be noted that, even then, as shown 
below, all static solutions to Eq. ([JJ can be obtained 
iteratively from Eq. |[7J) with an admissible K, starting 
from any admissible <po- 



C. Admissible values of integration constants 
and solutions from the nonlinear map 

Throughout this paper, without any loss of generality, 
we will always take A = ±1 and vary ft so as to vary A. 

Inadmissible values of integration constants, K and (f>o, 
are those for which, in Eq. JTJ), f(<j>o) < and 2 — A — 
A</>q = 0. For A = 1, we need to consider separately the 
regime of ft 2 < 2 (i.e. < A < 2) and the regime of 
very high discreteness, ft, 2 > 2 (i.e. A > 2). The third 
case is the case of A < where, however, the topology 
does not change as one goes from small to large h 2 . For 
A = 2 the two-point map Eq. ([7]) cannot be used in 
this form. This case is studied separately in Sec. Ill C 41 
It may be pointed out that similar studies have already 
been performed for the MC model in Q. 



1. Case of < A < 2 

The inadmissible regions of the plane (K, </>q) are 
shown by shaded areas in Fig. [JJ for different values of 
the discreteness parameter, indicated in each panel. The 
topology of the admissible regions is different for ft < 1 
and 1 < ft < y/2. We identify the portions of the ad- 
missible region occupied by different JEF solutions and 
mark the staggered solutions with " * " (see Sec. Ill Dp . 

Admissible regions for < A < 1 [see in (a) and (b) 
of Fig. [JJ are completely filled by the six JEF solutions: 
(sn/cn)*, (sndn/cn)*, sndn/cn, sn, 1/sn, and 1/cn (see 
Sec. HTDl) . 

In the range of 1 < A < 2 [see in Fig. [JJd)] the solution 
sndn/cn disappears, while the remaining five solutions 
persist and, there appears a portion of the admissible re- 
gion where we could not identify a JEF solution [marked 
with the question mark in (d)]. However, this solution 
can be easily obtained from the map Eq. JTJ and one 
finds that the solution is unbounded. 

The solution (sn/cn)* cannot be seen in (a)-(d) of Fig. 
[JJ because it occupies the range of — oo < K < K* with 
negative K* falling outside the frames of the panels. In 
(a) and (b) of Fig. [TJ the solution (sndn/cn)* takes a 
portion of K in between (sn/cn)* and sndn/cn solutions, 
but it falls out of the frames of these panels. 

The only bounded solution in this case is the sn solu- 
tion, which exists for < K < 1/2 when ft < 1 and in 
a narrower range of K otherwise. For example, in panel 
Fig. [TJd), the sn solution exists for 0.425 < K < 1/2. 
The kink solution, Eq. (|3"8]) , corresponding to K = 
and — 1 < 0o < 1, exists only for ft < 1 (i.e. only for 
< A < 1). 

Several static solutions generated by the map {7J for 
ft = 0.8 [see panel (b) of Fig. [JJ are presented in Fig. [21 
In Fig. HJa) we have the sndn/cn solution generated with 
K = —0.25, 0o = 0; in Fig. Hfb), the sn solution close to 
the hyperbolic limit generated with K = 10~ 8 , <fio = 0; 
in Fig. [DJe), the 1/sn solution close to the hyperbolic 
limit obtained for K = 10~ 8 , (f>o = 1.1; and in Fig. (2^d), 
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FIG. 1: Admissible and inadmissible (shaded) regions for the 
integration constants K and <j>o for A = 1 and different values 
of ft < \/2, indicated in each panel. The topology of the 
admissible regions is different for h < f and 1 < ft < y/2. 
We identify the portions of the admissible region occupied 
by different JEF solutions and mark the staggered solutions 
with " * " (see Sec. Ill D[) . The only bounded solution in this 
case is the sn solution which exists for < K < 1/2 when 
ft < 1 and in a narrower range of K otherwise. For example, 
in panel (d), the sn solution exists for 0.425 < K < 1/2. The 
kink solution, corresponding to K — and — 1 < cj>o < 1, 
exists only for ft < 1. The solution (sn/cn)* cannot be seen 
in (a)-(d) because it occupies the range of — oo < K < K with 
negative K being outside the frames of the panels. In (a) and 
(b), the solution (sndn/cn)* takes a portion of K in between 
(sn/cn)* and sndn/cn solutions, but it falls out of the frame 
of these panels. For the area marked with the question mark 
in (d) we could not identify a JEF solution. Several static 
solutions generated by the map of Eq. (J7| for ft = 0.8, panel 
(b), are presented in Fig. [2] 



the 1/cn solution constructed for K — 0.75, (f>o — 1.25. 
The only bounded solution is the sn solution shown in 
Fig. Hlb), when it assumes the form of alternating kinks 
and anti-kinks. For unbounded solutions, the amplitudes 
of some points are very large and they fall outside the 
frames of the figures. 



2. Case of A > 2 

The inadmissible regions are shown by the shaded ar- 
eas in Fig. [3] for different values of ft, indicated in each 
panel. These arc already quite different in structure 
from the inadmissible regions of Fig. [TJ notice in par- 
ticular the existence of such regions for positive K, near 
<f>Q = 0, while in Fig. [TJ admissible <po exist for any K. 
Notice also a qualitative change of the border between 
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FIG. 2: Static solutions generated by the map Eq. (0 for 
A = 1, ft = 0.8 [see panel (b) of Fig. [1] for different values of 
the integration constants, K and <f>o. We have: (a) sndn/cn 
solution (K = —0.25, <f>o =0); (b) sn solution close to the 
hyperbolic limit (K — 10 -8 , <f>o = 0); (c) 1/sn solution close 
to the hyperbolic limit (K — 10 -8 , 4>o = 1-1); and (d) 1/cn 
solution (K = 0.75, <j)o = 1.25). The only bounded solution 
is the sn solution shown in (b). Unbounded solutions have 
points with very large amplitudes falling outside the frames 
of the figures. 



admissible and inadmissible regions that takes place at 
h = 2 when the curvature of the border at the point 
(K, cj) ) = (1/2, 0) changes sign. One can see that, in this 
case, all static solutions are bounded. Admissible regions 
are completely filled by the three JEF solutions, cn, dn, 
and staggered dn, denoted as dn* (see Sec. Ill Dp . 

Examples of static solutions generated by the map of 
Eq. (0 for A = 1, h = 4 [corresponds to panel (d) of 
Fig. [3] for different values of the integration constants, 
K and cf>o, are plotted in Fig. 2] The solutions are: (a) 
dn solution (K = —0.05, (f>o — 1); (b) zigzag solution, 
which is actually the dn* solution in the limit of m = 
(K = 7/32, 0o = \/3/2); (c) dn* solution (K = 0.26, 
4>o = 1); and (d) cn solution with K{m) < (3h < 2K(m) 
close to the hyperbolic limit, where it assumes the form of 
the chain of staggered pulses, Eq. (O, (K = l/2+10~ 12 , 
0o = 1). 



3. Case of A < 

The inadmissible regions are shown by the shaded ar- 
eas in Fig. [5] for (a) h = 1.3 and (b) h — 4. In this case, 
the topology does not change as one goes from small to 
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FIG. 3: Same as in Fig. [JJbut for a regime of high discrete- 
ness, h 2 > 2. A qualitative change of the border between ad- 
missible and inadmissible regions takes place at ft = 2 when 
the curvature of the border at the point (K,(f>o) = (1/2,0) 
changes sign. All static solutions are bounded. Admissible 
regions are completely filled by the three JEF solutions, cn, 
dn, and staggered dn, denoted as dn* (see Sec. IIIDfl . Several 
static solutions generated by the map of Eq. (0 for h — 4, 
panel (d), are presented in Fig. [4] 



large h 2 . All static solutions are bounded. The admissi- 
ble regions are completely filled by the three JEF solu- 
tions, cn, dn, and staggered dn, marked as dn* (sec Sec. 

EH. 

Several static solutions generated by the map of Eq. 
for h = 1.3 [refer to panel (a) of Fig. [5] are presented 
in Fig. El In (a), K = -0.5 - 1CT 8 , = 1; in (b), 
K = -0.5 + 10~ 8 , </-o = 1; in (c), K = -0.2, </> = 1; and 
in (d), K = 5.3, <j)o — 2. The solutions in (a) and (b) are 
the cn and dn solutions close to the hyperbolic function 
limit, respectively. These solutions can be regarded as a 
lattice of pulses, Eq. (|40|) . Solutions shown in (c) and 
(d) are dn and staggered dn solutions, respectively. 



4. Case A = 2 

In this case, the two-point map of Eq. (O reduces to 
€ (€€+1 - 2^>„0„+i +K)=0, (16) 

from which one has the trivial solution <p n = and, for 
<j) n 7^ 0, one finds 



n+1 = -L[K± X /K(K- I 



(17) 



This map defines the two-parameter space of solutions. 
Admissible regions of the (K,<j>o) plane correspond to 
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FIG. 4: Static solutions generated by the map of Eq. Q for 
A = 1, h = 4 [see panel (d) of Fig. [3] for different values 
of the integration constants, K and <j>a. All solutions are 
bounded. We have: (a) dn solution (K — —0.05, <f>a = 1); 
(b) zigzag solution, or the dn* solution in the limit of m = 
(K = 7/32, 4> = v/3/2); (c) dn* solution (K = 0.26, <f> = 1); 
and (d) the cn solution with K(m) < f3h < 2K(m) close to 
the hyperbolic limit, where it obtains the form of the chain of 
staggered pulses, Eq. (gSJ, (K = 1/2 + 10~ 12 , <po = 1). 



K < and K > 1. Due to the symmetry of Eq. (|16p . 
one can interchange <p n and in Eq. (|17p . From 

the two different roots of Eq. (fTT]) one must take 4>n+i 
such that the three point problem Eq. ([T]) is satisfied for 
4> n ; this would guarantee that one keeps moving along 
the same solution branch. It is sufficient to take 4> n +i 
different from <f> n -i, if the roots are different. 

It is convenient to take the integration constant in the 
form 



K 



(1 



4a 



(18) 



and consider — 1 < a < 1 with a / 0. With this choice, 
for —1 < a < we have — oo < K < 0, while for < a < 
1 we have 1 < K < +oo. 

Taking for the initial value (f>o = b, where b =/= is 
positive or negative real number, one can obtain from 
the map of Eq. (fTT)) the following sequence, that satisfies 
both Eq. flTJ and Eq. (JTBJ), 



(1 + a) b (l + a)a b (l + a)a 2 



26 



26 



26 



.(19) 



Note that for one of the bounds of the existence region, 
K = 1 (a = 1), solution Eq. (fT9")l gives the following 
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FIG. 5: Same as in Fig. [Qand Fig. |3]but for A = -1. In this 
case, the topology does not change as one goes from small to 
large h 2 . All static solutions are bounded. The admissible 
regions are completely filled by the three JEF solutions, cn, 
dn, and staggered dn (marked as dn*). Several static solutions 
generated by the map of Eq. (0 for h = 1.3, panel (a), are 
presented in Fig. [6] 
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family of two-site periodic solution 



(20) 



On the other hand, the other bound of K = (a = 
— 1) gives from Eq. (fl9|) the following four-site periodic 
solution 



b„ : (...,6,0,-6,0,6,0,...). 



(21) 



It is worth pointing out that for the MC model (U]) at 
A = 2, the two-point map has the form 

(i - ^M+i(i - <t>l) - (i - €) + c] = o, (22) 

and, for (j) n ^ ±1, we obtain 



= ±\ 1 



c 



1 



(23) 



Based on this expression, starting from any admissible 
value of (f>o and for any chosen C, one can construct a so- 
lution iteratively. The topology of the admissible regions 
in this case is not as simple as in the case of model Eq. 
(□J at A = 2. 

A four-site periodic solution for <f) n in this case is given 

by 



4> n : (..., a, 6, —a, —6, a, 6, —a, —6, a, ...) 
C=(l-a 2 )(l-6 2 ), 



(24) 



where a, 6 are any real numbers, positive or negative. In 
the special case of C = 0, apart from the above solution, 
another solution for the sequence <p n is the arbitrary se- 
quence of ±1 with put at random between —1 and +1 
or between +1 and —1. 



FIG. 6: Static solutions generated by the map of Eq. (0 
for A = —1, h — 1.3 [see panel (a) of Fig. [5] for different 
values of the integration constants, K and <j>o. In (a), K = 
-0.5 - 10" 8 , (j> = 1; in (b), K = -0.5 + 10" 8 , <j) = 1; in 
(c), K = -0.2, 4>a = 1; and in (d), K = 5.3, 4> = 2. All 
the solutions are bounded. The solutions in (a) and (b) are 
the cn and dn solutions close to the hyperbolic function limit, 
respectively. The solution in (c) is the dn, while in (d) we 
show the staggered dn solution. 



D. Static Jacobi elliptic function solutions 

Motivated by the results of the previous section, let us 
now proceed to analytically identify several static JEF 
solutions to Eq. (JTJ) . The solutions can be derived from 
the three-point static problem of Eq. ([T]) with 4> n = 0. 
In this case the use of the JEF identities listed in the 
Appendix II is particularly helpful. Alternatively, solu- 
tions can be derived from the DFI of Eq. ([7]). In this 
case, one also obtains the relation between the integra- 
tion constant K and the parameters of the JEF solution 
thus establishing a link between the solutions obtained 
by two different approaches. 

The results will be presented in the following way. We 
will first describe the general form of the derived JEF so- 
lutions, then give the parameters for the particular ones, 
and finally, we will characterize the solutions and describe 
the ranges of parameters where different solutions exist. 

The general form of the solutions is 

4> n = ±S , Asn p (Z,m)cn 9 (Z,TO)dn r (Z,m), 
Z = (3h(n + x ), (25) 

where sn, cn, and dn are the Jacobi elliptic functions, 
< m < 1 is the modulus of JEF, A and /3 are the 
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parameters of the solution, and xq is the arbitrary initial 
position. For non-staggered solutions S = 1, while for 
staggered ones S = (— 1)™. Finally the integers p,q,r 
specify a particular form of the solution. 

For the sake of brevity we introduce the following no- 
tations 

s = sn(/3/i, m), c= cn(f3h,m), d = dn{j3h,m). (26) 

Particular solutions have the following form and are 
characterized by the following parameters: 
The sn solution, (p, q, r) = (1, 0, 0), 



A(l - ms 4 ) 



.o (2 — A)ms 2 , 2ms cd 
A 2 = ± ^ , A - 

K = —, 5=1. 

m 



The cn solution, (p, q, r) = (0, 1, 0), 

2 (A — 2)ms 2 .n 2ms 2 c 



A 2 = 



Ad 2 



A = — 



A 4 



K = -(l-m) — , 5 = 1. 
m 

The dn solution, (p, q, r) = (0, 0, 1), 



.4 



(A-2)s 2 



A 1 



2s 2 d 



A2= (2-A)ms 2 a2 



A 



A 2 = 



A(l - ms 4 ) 



.4 



K=— , 5=(-l)". 
m 

The staggered cn* solution, (p, q, r) — (0, 1, 0), 



A 



(A - 2)ms 2 



Ad 2 



A 



2ms 2 c 



A{d 2 — ms 2 c 2 ) 



A 4 



K = -{l-m) — , 5=(-l) n . 
m 

The staggered dn* solution, {p, q, r) = (0, 0, 1), 



A 



2 _{A-2)s 2 



A 2 



2s 2 d 



A 2 — a 2 

A(l - ms 4 ' A 

if = ^ 4 m, 5=1. 



(27) 



A(d 2 — ms 2 c 2 ) 



(28) 



Ac 2 ' Ji A(c 2 -s 2 d 2 )' 
iY=(l-m)A 4 , 5 = 1. (29) 

The staggered sn* solution, (p,q,r) — (1,0,0), 

2ms 2 cd 



(30) 



(31) 



Ac 2 ' c 2 — s 2 d 2 ' 

K = (l-m)A 4 , 5=(-l) n . (32) 

The 1/sn solution, (p,q 7 r) = (—1,0,0), 

2s 2 cd j9 (2-A)s 2 



(33) 



The 1/cn solution, (p,q,r) — (0,-1,0), 



A = 



2(1 - m)cs 2 
A{d 2 — ms 2 c 2 



A 2 = 



(2- A)(l- m)s 2 
Ad 2 



.. 5=1. 



(34) 



1 — m 

The staggered (sn/cn)* solution, (p,q,r) = (1,-1,0), 



A = 



K 



(l-m)(2-A)f 



.4 4 



Ac 2 

, 5 = (-iy 



A = - 



2(1 - m)s 2 d 
A(c 2 - s 2 d 2 ) : 



1 — m 

The sndn/cn solution, (p,q,r) = (1,-1,1), 



(35) 



a2 _ 2s 2 d 2 {l- 2ms 2 + ms 4 ) ^ 2 _ (2 - A)s 2 d 2 



A(c 4 - s 4 d 4 ) 



Ac 2 



K = A\ 5 = 1. 



(36) 



The staggered (sndn/cn)* solution, (p,q,r) 
(1,-1,1), 



2s 2 d 2 (l - 2ms 2 + ms 4 ) ^ _ (2 - A)s 2 rf 2 
A(c 4 - s 4 d 4 ) 



A 2 = 

K = A 4 , S = (-1)". 



, A 2 



Ac 2 



(37) 



It is worth making the following remarks: 
To specify a particular solution, for given model pa- 
rameters h 2 and A, and for chosen m, one has to solve 
the first two equations entering Eqs. (f2"7) - (f3~T|) in or- 
der to find the parameter f3 and then the amplitude A. 
If the equation for (3 has no solutions in the range of 
< /3h < 2K(m), this implies that for chosen h 2 , A, and 
m the corresponding JEF solution does not exist. Some 
of the roots for (3 can give imaginary amplitude A and we 
reject such solutions. On the other hand, if one or more 
roots for (3 in the range of < (3h < 2K{m) are found, 
then a JEF solution can be constructed for each root 
provided that the corresponding amplitude A is real. In 
TableU for different JEF solutions, we summarize the lo- 
cation of roots for the parameter /3, corresponding to real 
amplitude A and give the range of parameter A where 
the solutions exist. For comparison, similar results for 
bounded JEF solutions to the MC model of Eq. (TJ| are 
given in Table ITT1 

In each of Eqs. (f!Z7| - (f3"T]) we give the corresponding 
expressions for K in terms of the parameters of JEF so- 
lutions. These expressions link the JEF solutions to the 
nonlinear map Eq. ((?]), so that for each JEF solution one 
can find the corresponding portion of the (K,<f>o) plane. 
For example, the sn solution exists for < A < 2 (see 
Table [TJ) and this case is shown in Fig. [T] Then, the 
sn solution having modulus m and amplitude A, in the 
(K, 4>o) plane will have abscissa K that can be found for- 
given K from Eq. © and the ordinate will vary in the 
range — A < <fio < A. Considering different values of 
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TABLE I: Range of parameter A and location of the root 
for parameter j3 corresponding to real amplitude A for JEF 
solutions to Eq. {1} 



Solution 


Root for /3 within 
< ph < K(m) 


Root for P within 
K(m) < ph < 2K(m) 


sn, Eq. J27J) 


< A < 2 


no solution 


cn, Eq. Q28J 


A < 


A > 2 


dn, Eq. (|29J) 


A < 


A > 2 


* J!Wt 1 

sn', Eq. Q30J) 


no solution 


< A < 2 


cn*, Eq. (|3T| 


A > 2 


A < 


dn*, Eq. (|32J 


A > 2 


A < 


1/sn, Eq. ([33| 


< A < 2 


no solution 


1/cn, Eq. (|34|l 


< A < 2 


no solution 


(sn/cn)*, 
Eq. (ESJ 


< A < 2 


< A < 2 


sndn/cn, 
Eq. fl36]) 


< A < 1 


< A < 1 


(sndn/cn)*, 
Eq. §7} 


< A < 2 


< A < 2 



TABLE II: Same as in Table U but for the bounded JEF so- 
lutions to the MC model Eq. @ 



Solution 


Root for /3 within 


Root for f3 within 




< (3h< K(m) 


K(m) </3h< 2K(m) 


sn 


< A < 2 


A > 2 


cn 


A < 


no solution 


dn 


A < 


no solution 


sn* 


A > 2 


< A < 2 


cn* 


no solution 


A < 


dn* 


no solution 


A < 



< m < 1, one can see that the sn solution occupies the 
portions marked with "sn" in Fig. [T] 

The solutions presented by Eqs. (f2"T)) - (|32p are bounded 
while the other solutions are unbounded. 

The sn solution and its staggered form sn* are in fact 
equivalent in the sense that any sequence of </>„ produced 
by one of them can also be produced by the other. As 
a result, sn and sn* solutions occupy the same portion 
of the (K, (f>o) plane (see Fig. QJ. The same can be said 
about the pair of solutions cn and cn*. However, the dn 
solution is essentially different from its staggered form, 
dn*, and also the sndn/cn solution is essentially different 
from its staggered form, (sndn/cn)*. 

E. Hyperbolic function static solutions 

In the limit of m — > 1 the sn solution given by Eqs. 
(f25|) and (|27|) reduces to the kink (antikink) solution 

<p n = ± tanh\ph(n + x )], (38) 

where xq is an arbitrary shift, (3 satisfies 

= (39) 

and K = l, i.e., K = 0. 



TABLE III: Range of parameter A for hyperbolic solutions to 
Eq. © 



Solution 


Exists for 


Kink, Eq. O 


< A < 1 


Pulse, Eq. (140 j| 


A < 


Staggered kink 


no solution 


Staggered pulse, Eq. (1421) 


A > 4 



TABLE IV: Same as in Table [TTT] but for MC model Eq. © 



Solution 


Exists for 


Kink 


< A < 2 


Pulse 


A < 


Staggered kink 


A > 2 


Staggered pulse 


no solution 



In the limit of m — > 1, both the cn solution and the 
dn solution with parameters given by Eq. (|28[) and Eq. 
(f2"9")l . respectively, reduce to the same pulse solution 

4> n = ±Asech[/3h(n + x )], (40) 

where Xq is an arbitrary shift, (3 and A can be found from 

2 - A 

cosh(/%) = — ^— , 

KA 2 = (A-2)sinh 2 (/%), (41) 

and K = 0, i.e., for A = -1, K = -1/2. 

We now turn to the case of staggered hyperbolic func- 
tion solutions. It is easily checked from Eq. ([3TJ]) that 
the TI model of Eq. {1} does not support the staggered 
kink solution (even though it supports the sn* solution). 
However, it supports the staggered pulse solution 

4> n = ±(-l)Msech[/3/i(n + x )], (42) 

where Xq is an arbitrary shift, j3 and A can be found from 

A - 2 

cosh(/3/i) = — — , 

KA 2 = (A-2)sinh 2 (/3fr). (43) 

and K = 0, i.e., for A = 1, K = 1/2. 

In Table HTT1 we present the ranges of model parameter 
A where the hyperbolic function solutions exist in the 
model of Eq. {l} and, for comparison, in Table IIVI wc 
give similar results for the MC model of Eq. |4]) . 

III. LINEAR STABILITY AND DYNAMICS 
FOR BOUNDED SOLUTIONS FOR h 2 < 2 

A. Linear stability 

Introducing the ansatz (f> n (t) — + e n (t) (where 0° 
is an equilibrium solution and s n (t) is a small perturba- 
tion), we linearize Eq. (Q]) with respect to e n , in order 
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to examine the linear stability of the above-established 
solutions. In this way, we obtain the following equation: 



Ae„ 



-A<^°+lEn-l ~ A0°_ 1 0° + 1 e„ - A$_i$Un+i- (44) 

For the small-amplitude phonons, s n — exp(ikn + iuit), 
with frequency to and wave number fc, Eq. (|44|) is reduced 
to the following dispersion relation: 



+A0° ( />° +1 e- jfc + A^_ 1 i 



aO „ifc 



(45) 



where it is assumed that 0° is a uniform steady state solu- 
tion (or a zigzag solution; see below). The corresponding 
eigenvalue problem has a non-symmetric matrix. This 
fact is a manifestation of the non-Hamiltonian nature of 
Eq. 0. 

Stability of vacuum solutions. From Eq. (|45[) . the spec- 
trum of the vacuum solution 0" = ±1 (at A = 1) is 



uo z = 2 



4 , 9, 9 /fc 



(46) 



This solution is stable for h 2 < 2. 

The spectrum of the vacuum solution 0° = (at A 
-1) is 



w =1 + 7? sin 2 



(47) 



and this solution is stable for any value of h 2 . 

Stability of kink solution. As it can be seen from Eq. 
(|39f , the kink solution Eq. (|38|) exists only for A = 1 and 
h 2 < 1 . The spectrum of lattice with a kink as a function 
of lattice spacing h is shown in Fig. [7] for the on-site 
kink, x — (open circles), and inter-site kink, xq = h/2 
(dots). The corresponding kink profiles are shown in the 
inset for h — 0.8. Two solid lines show the borders of 
the spectrum of vacuum, Eq. (l4"6")) . The kink possesses 
the translational internal mode with ui 2 = and also the 
localized modes with frequencies below the phonon band 
and a localized mode with frequency above the phonon 
band. The kink is stable in the whole range of existence, 
h 2 < 1. Notice also that as the continuum limit h — ► 
is approached, a number of internal modes that exist for 
larger h disappear into the phonon band of Eq. (1461) ; the 
only internal mode that persists with oj — > \/3, as h — > 
is the well-known continuum </> 4 model internal mode (see 
e.g. ((J and references therein). 

Instability of the zigzag solution. In the sinusoidal func- 
tion limit (m — 0), the staggered dn solution defined by 
Eq. (|25|) and Eq. (|32jl reduces to the zigzag solution 



4> n = (-i)M, 



A = 



A-4 
A ' 



(48) 



Note that this solution is valid if A > 4 or if A < 0. 
Substituting this solution into Eq. (|4"5")) we obtain the 



dispersion relation u?h 2 = 4 - 2 A + 4(A - 3) sin 2 (£;/2). 
One can see that the spectrum of the zigzag solution 
always contains imaginary frequencies in both cases, A > 
4 and A < 0, and thus, the solution is always unstable. 

Stability of solutions at A = 2. In this case (see Sec. 
Ill C 4|) . we have two bounded solutions, Eq. (f20|) and 
Eq. ([2l"|). For the two-periodic solution Eq. (|2"0)). the 
vibrational spectrum has two branches, (/i 2 w 2 )i,2 = a 2 + 
1/a 2 ± v /(a 2 + l/a 2 ) 2 -2 + 2cosfc. Since a 2 + l/a 2 > 2 
for any a, uj 2 is always non-negative and the solution 
is always stable. Numerically we found that the four- 
periodic solution Eq. (|2"Tj) is always unstable. This can 
be justified by introducing a small deviation e into the 
position of a particle with <p = 0, while keeping all other 
particles in their equilibrium positions. The right-hand 
side of Eq. ([XJ) gives the force acting on the particle equal 
to 2ea 2 . Since the sign of the force coincides with the sign 
of the displacement, the solution is indeed unstable with 
respect to the considered perturbation. 

For the MC model Eq. (UJ) in the case of A = 
2, we found that the solutions ...,—1,1,-1,1,... and 

— 1, 0, 1, — 1, 0, 1, ... are unstable. This can be checked 
by introducing a small displacement into position of a 
particle with <f> — 1, keeping other particles at their equi- 
librium positions. Wc then calculate the increment of 
force acting on this particle and find that it has the same 
sign as the displacement, which proves the instability 
with respect to this type of perturbation. 

Calculating the increment of force acting on a particle 
with unperturbed amplitude <f> = 1 in response to a small 
displacement e of this particle, we found that it has the 
same sign as the displacement. Thus, these two solutions 
are unstable and this fact was also confirmed numerically. 

Stability of other solutions. We also examined numer- 
ically the stability of other ones among the above-found 
solutions. We studied bounded, periodic Jacobi elliptic 
function solutions depicted in Fig. [T]and Fig. EJa), i.e., 
the sn, cn, dn, and staggered dn solutions. The length 
of the chain was compatible with the period of the so- 
lution and it contained a few periods. The eigenvalue 
problem obtained by replacing the left-hand side of Eq. 
(|44p with — u> 2 e n was solved numerically under periodic 
boundary conditions. The sn, cn, dn, staggered dn, and 
the pulse solutions were found to be generically unstable, 
i.e., having negative eigenfrequencies a; 2 in their numer- 
ical spectra. 

It is important to note that all the solutions possess 
the zero-frequency translational eigenmode due to the 
TI nature of the lattice of Eq. |T]). As an overview of the 
results of the stability analysis, we find that, among the 
bounded solutions, only the kink solution, the vacuum 
solutions, and the solution of Eq. (T5U|) (A = 2 case) can 
be stable. Note that exactly the same conclusion was also 
reached in the case of the MC model [9j in case A ^ 2. 
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FIG. 7: Frequencies of the kink's internal modes for different 
magnitudes of the discreteness parameter h. Results for the 
on-site (circles) and inter-site (dots) kinks. The corresponding 
kink profiles are shown in the inset for h = 0.8. Two solid 
lines show the borders of the spectrum of vacuum, Eq. (|46[) . 
The kink possesses the translational internal mode with ui 2 = 
and also the localized modes with frequencies below the 
phonon band and a localized mode with the frequency above 
the phonon band. The kink is stable in the whole range of its 
existence, h 2 < 1. 



B. Dynamics of unstable solutions 

In the figures I5HTTH we give examples of the dynamical 
behavior of the unstable static solutions perturbed at t = 
by adding a small displacement to the 0-th particle. 

In Fig. [5] the sn solution was perturbed at t = by 
shifting the position of n = particle by — 10~ 2 . The 
model parameters are chosen as A = 1, h = 0.25 and the 
parameters of the sn solution are m ~ 0.777, xo = 0. 
Up to the time t w 45 the long-wave mode of devia- 
tion from equilibrium, growing exponentially with time, 
remains indistinguishable in the scale of the figure. Af- 
ter that, the sn solution gradually slides down into the 
potential well at 4> = — 1 and a set of breather- like os- 
cillating pulses is formed. The oscillations are accompa- 
nied by phonon radiation and their amplitude decreases 
in time. The instability picture presented in Fig. [8] can 
also be interpreted in terms of the kink-antikink interac- 
tions. Such interactions are attractive and the chain of 
mutually attractive quasi-particles, as it is well-known, is 
unstable and cannot remain equally spaced in the pres- 
ence of small perturbations. 

Similar behavior is demonstrated by the unstable pulse 
perturbed at t — by shifting the position of the n = 
particle by — 10~ 4 (see Fig. The pulse slides down 
into the potential well at (f> — and oscillates by radiating 
energy. The pulse is shown at the moment of maximum 
displacement of the central particle. The period of os- 
cillation of the pulse is roughly equal to 7. The model 



FIG. 8: Time evolution of the sn solution perturbed at t = 
by shifting the position of n — particle by — 10 -2 . The 
model parameters are F = A = 1, h = 0.25 and the parame- 
ters of the sn solution are m ~ 0.777, xo = 0. 



parameters are A = —1, h = 0.25 and the pulse is placed 
at xq = 0. We note in passing that the instability of 
the single pulse has been demonstrated quite generally 
for discrete TI settings in the case of the (f> 4 equation in 
@ ; on the basis of that proof, we expect the individual 
pulses (as well as solutions consisting of concatenations 
of such pulses) to be unstable in the present model as 
well. 

In Fig. [10] the time evolution of the unstable dn so- 
lution perturbed at t — by shifting the position of the 
n = particle by — 10 -4 is presented. The model param- 
eters are A = —1, h = 0.25 and the parameters of the 
dn solution are m — 0.3, xo = 0. The dn solution has 
particles situated near the top of the potential barrier 
at <fi = 1 and, being perturbed as described above, the 
particles, one after another, slide down into the potential 
well at 4> = and oscillate near the potential well. 

Note that the perturbation of the opposite sign (out- 
ward of the potential well at c/> = 0) would cause the 
divergence of the pulse and the dn solutions shown in 
the top panels of Fig. [5] and Fig. QJJ1 respectively. We 
note that qualitatively similar behavior was also observed 
in the case of the MC model [9(. 



IV. DISCUSSION AND CONCLUSIONS 

For a non-standard, translationally invariant discrete 
lattice, namely the discretization of the <j) A model of Eq. 
(UJ, we have reported a number of important develop- 
ments. We have explicitly identified the momentum con- 
servation law given by Eq. ([5]) , produced a set of static 
solutions in terms of the Jacobi elliptic functions, includ- 
ing the staggered ones, see Sec. HID I and their hyperbolic 
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FIG. 9: Time evolution of the pulse solution perturbed at 
t = by shifting the position of n — particle by — 10~ 4 . 
The pulse is shown at the moment of maximum displacement 
of the central particle. The period of oscillation of the pulse is 
roughly equal to 7. The model parameters are F = A = — 1, 
h — 0.25 and the pulse is placed at xq = 0. 



of solutions of the static problem, and identify the asso- 
ciated conservation law for such discrete, translationally 
invariant models. 

An important feature of our approach is that from the 
nonlinear map Eq. ([7]), any static solution of Eq. (fTJ) 
can be obtained recurrently, solving at each step a sim- 
ple algebraic problem, in our case, a quadratic equation. 
Admissible values of the integration constants K and 0o 
have been plotted for positive and negative A and for dif- 
ferent values of the discreteness parameter h 2 , see Figs. 
[1] [3l and \5[ Any point inside the admissible region cor- 
responds to the static solutions of Eq. (JTJ) . 

We expect that the results presented herein for the 
discrete (j) 4 equation can be applied to analogous discrete 
nonlinear Schrodinger type equations, similarly to what 
was discussed in the recent works 112, [11 SSI. An- 
other very natural extension of the present approach, de- 
tailing how to obtain an analytical handle on solutions of 
such translationally invariant models, would be to con- 
sider higher dimensional settings. Very little is known 
even in 2+1-dimensional spatially discrete models about 
analytical solutions, hence such investigations would be 
a particularly interesting next step in this direction of 
research. 
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FIG. 10: Time evolution of the dn solution perturbed at t = 
by shifting the position of the n = particle by — 10~ 4 . 
The model parameters are F — A = — 1, h = 0.25 and the 
parameters of the dn solution are m = 0.3, xq — 0. 



limit. We have fully analysed the existence of such so- 
lutions, as well as detailed their linear stability. In Ap- 
pendix I, we have demonstrated how, using a JEF solu- 
tion of Eq. (fT]), one can construct the nonlinear map (or 
DFI) of Eq. ||7J). Importantly, we have proposed this as 
an illustrative example of how to formulate the discrete 
first integral approach, how to link it to the Jacobi elliptic 
function solutions, how to construct the complete space 



Appendix I 

In this Appendix we show how the two-point maps can 
be constructed directly from the known exact solutions. 
For example, for the tanh solution Eq. one can find 
4> n -i or 4>n+i using the identity tanh(a: + y) — (tanhx + 
tanhy)/(l + tanh x tanh y) and then, making use of Eq. 
(155)) and Eq. write down the two-point map that 

generates the tanh solution 



£>n±l 



4>n±T 

l±T<j> n 



T = 



A 



2- A 



(49) 



This map does not involve the integration constant be- 
cause for the tanh solution we have already set m = 1. 
The initial value can be chosen arbitrarily from the range 
-1< 00 < I- 

Similarly, for the sn solution, Eq. (|25|) and Eq. (f27|) . 
using the formula for sn(x + y), we can write 



.4 



cdS ± sCD 
1 — ms 2 S 2 ' 



(50) 



where we have denoted 



S = sa.\(3h{n + xq), m], 
C = cn\j3h{n + xo), m], 
D = dn\flh(n + xq), to] 



(51) 



and also used the notations of Eq. (|2"rj)) . According to 
Eq. ([2T)l we have cj> n = AS and, using the identity for the 
Jacobi elliptic functions C 2 D 2 = I — (I + m)S 2 + mS 4 , 
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we rewrite Eq. (|50[) in the form of the two-point map 



^" ±1 = A-msWJA ■ (52) 

The integration constant here is < m < 1. As soon as 
to is specified, one can find (3 and A from Eq. ([27]) and 
then s, c, d from Eq. (|26p . thus defining the map Eq. 
(|52p . The corresponding sn solution can be constructed 
by iterating Eq. ((52)) starting from any initial value taken 
from -A < 4>o < A. 

By analogy, two-point maps can be constructed for any 
solution given in Sec. Ill Dl and Sec. Ill El However, each 
of the maps constructed in this way is parameterized by 
the modulus to and can be applied only to the corre- 
sponding solution. In our case, it is possible to substi- 
tute to with a universal integration constant so that the 
resulting two-point map can generate all static solutions 
of Eq. (TT]). The universal integration constant can be 
taken, e.g., as 



2 V m J 



then we have 

k = i 

From Eq. (|2"T|) one can derive 



2K 
X 



A^_ 

to 



cd 



(1+m) A 2 
Z, ± '- = A, 



m 



(53) 



(54) 



(55) 



where Z and X are given by Eq. (fT0|) and Eq. (TTTj) , re- 
spectively. With the help of Eq. (27]) and Eqs. (I53l- (j55l) 
we transform Eq. (fB"2"|) to the form of Eq. ([7]) which does 
not include explicitly the parameters of the sn solution. 



Appendix II 



TO,cn(a; + a)cn(x)cn(x — a) 
= — ds(a)ds(2a)[cn(x + a) + cn(x — a)] 
+ds 2 (a)cn(x) , 



(57) 



TOdn(x + a)dn(x)dn(a; — a) 

= — cs(a)cs(2a)[dn(x + a) + dn(x — a)] 

+cs 2 (a)dn(x). (58) 
On the other hand, in the case of the MC model 0], 
the following identities are required for deriving the JEF 
solutions Hal 



TO.sn 2 (a;)[sn(a: + a) + sn(x — a)] 
= ns 2 (a)[sn(a; + a) + sn(x — a)] 

— 2cs(a)ds(a)sn(x), (59) 



TO.cn 2 (a;)[cn(a: + a) + cn(ir — a)] 
= — ds 2 (a)[cn(x + a) + cn(x — a)] 

+2cs(a)ns(a)cn(a;), (60) 



TOdn 2 (a;)[dn(a; + a) + dn(x — a)] 
= — cs 2 (a)[dn(a; + a) + dn(x — a)] 

+2ds(a)ns(a)dn(a;). (61) 



In the derivation of the Jacobi elliptic function solu- 
tions presented in Sec. Ill Dl the following identities for 
Jacobi elliptic functions have been used [25| 
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